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• 1 INTRODUCTION 

o 

^ I ' The existence of a supermassive BH in the centers of 
I I galaxies is well-established (iKormendv fc Richston3 1 19951 : 
. . i lRichstone et al.lll998l : lKormendv fc Gebhardtll200ll '). Strong 
^ ■ links exist between the supermassive BH and host galaxy 
properties, as evidenced by the BIf mass-stellar velocity 
. dispersion (Mbh — g ) relation (|Ferrarese fc MerrittI l2000l : 
■ iGebhardt et al.llioool '). amongst others. The Mbh — cr rela- 
tion has been demon strated not to be a selection effect by 
iGiiltekin et al.l (|201ll '). to be the strongest co rrelat i on be - 
tween Mbh and galaxy properties by Beifioii et al.l (l201lf) 
and h as been estimated most recently by IGiiltekin et al.l 
(|2009t ) based on 49 Mbh measure ments and 19 Mbh up- 
per limits, bv lGraham et al" ||2011) bas ed on 64 Mbh mea- 
surements, and bv lBeifiori et al.l ( 201ll ) based on 49 Mbh 
measurements and 94 Mbh upper limits. 

Despite the wealth of observational data, the origin 
of this relation is not firmly established. There are several 
theoretical models that explain the origin of the observed 
correlations between BIf mass and galaxy properties. 
Feedback from BH accretion on the hosting galax:y is one 



ABSTRACT 

We investigate the evolution of the Mbh — ct relation by examining the relationship 
between the intrinsic scatter in the A/bh — cr relation and galaxy bolometric nuclear 
luminosity, the latter being a probe of the accretion rate of the black hole (BH). Our 
sample is composed of galaxies with classical bulges when possible, of which 38 have 
dynamically measured BHs masses, and 17 have BHs masses measured by reverbera- 
tion mapping. In order to obtain the bolometric nuclear luminosity for galaxies with 
low nuclear luminosity, we convert the X-ray nuclear luminosity measured by Chandra 
to bolometric luminosity. We find that the scatter in the Mbh — c relation is uncorre- 
lated with nuclear luminosity over seven orders of magnitude in luminosity, with the 
high luminosity end approaching the Eddington luminosity, LEdd- This suggests 
that at the present epoch galaxies evolve along the A/bh ^ c relation. This conclusion 
is consistent with the standard paradigm that BHs grow contemporaneously with their 
host stellar spheroids. 



accretion - galaxies: nuclei - galaxies: kinematics 



further fueling of the BH i Kauffmann fc Haehnelt' '200g|: 
Wvithe fc Loebll2003l: iDi Matt co ct al. 2005; Murray ct ^ 
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Hopkins et al 



Robertson et al 
20091 ; 



Ijohansson et al.l l2009l : 
Silk fc Nusseil I2OI0I '). This feed- 
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prop osal (|Silk fc Rees .199^: iFabianI Il999l : iBurkert fc "siilj 
I2OOII ). Simulations often involve galaxy mergers with strong 
inflow of gas that feeds the BH, powers the quasar and 
expels enough gas to quench both star formation and 



back from active galactic nuclei (AGN) regulates the 
BH-galaxy systems, and leads to tight BH mass-galaxy 
property relations. This scenario predicts that Mbh and a 
coevolve. 

Other mechanisms have been proposed to explain the 
tight corre l ation betw een BH mass and galax y proper- 
ties. IPend (|2007l 'l and iHirschmann et all (|2010l ') construct 
a model for the origin of the Mbh — Mi,uige relation in 
which mergers do not lead to accretion-based growth on 
the BH. In this model a tight A/bh — Mbuige relation is 
established throu g h the central limit theorem. Recently, 
IJahnke fc Macciol (|201ll ) build on this model by including 
BH accretion and star formation (these processes are un- 
correlated in the model). They conclude that a causal link 
between galaxy growth and BH growth is not necessary for 
obtaining the observed BH mass-galaxy property relations. 

In this article we investigate how galaxies evolve in the 
Mbh — o" plane, and thereby place constraints on these and 
other models for the origin and evolution of the Mbh — o" 
relation. We focus on local galaxies with classical bulges, 
and investigate the scatter in the Mbh — cr relation as a 
function of galaxy nuclear luminosity. Our results indicate 
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a scenario where BH mass and a evolve along the Mbh — cr 
relation, thereby favoring models where BH mass and host 
stellar spheroids coevolve. 



2 METHOD 

First, we describe the criteria adopted to select our sam- 
ple and the methods used to analyze the data. We con- 
sider the samples co mpiled by iGiiltekin et al.l l|2009l ) and 
iGraham etlo] l|201ll ). which include galaxies with dynami- 
cally measured BHs and stellar veloci ty dispersions . In ad - 
dition, we include galaxies studied bv I Greene et all ||2010D . 
who measured the central BH mass via masers. We refer to 
the sum of these samples as our dynamically-based sample. 
We consider a second set^f galaxies with BHs estimated via 
reverberation mapping (jWoo et al,.2010, ). 

We select galaxies w ith classical bulg e s wh en possi- 
ble. This is done because iKormendv et all |20l3) observe 
that Mbh does not correlat e with the p roper ties of galaxy 
disks or pseudobulges, and ISani et al.1 (j201lf ) find smaller 
intrinsic scatter of BH mass-host galaxy property relations 
when excluding galaxies with pseudobulges. According to 
bulge classification of the galaxies included in iGreene et al.l 
l|2010t ). whose classification relies on the galaxy morphology 
and stellar population property, we select the only galaxy 
(N1194) that has a classical bulge with nuclear luminosity 
at lO"^'^'^ LEdd, and we also select the galaxy (UGC 3789), 
which has an unclassified bulge with nuclear luminosity at 
-|^g-o.82 j^^^^ poj- the other dynamically- based galaxies, w e 
select classical bulge galaxies according to (|Sani et al.ll201l[ ) , 
who classify galaxies with classical bulges by selecting galax- 
ies which have Sersic indices higher than two. We include all 
the reverberation-based galaxies as it is difficult to classify 
the morphology of galaxies with AGN. 

Next, we estimate the b olometric nuclear lu minosity. 
First we consider the sample o fin reene et al.l(|2010l ). who es- 
timate the nuclear bolometric luminosity for their sample us- 
ing O [HI] , which is strong a nd ubiquitous in obscured mega- 
mase r systems ( Kauffmann et al.l l2003l : IZakamska et al.l 
I2OO3I ). In this approach, 0[III] luminosity is converted 
to M2500, where M2500 is the magnitude at 2500 A, 
and the n to bolometric l uminosity following the M2500 — 
LfOIII] jReves et a l. 20081 and M2500 bolometric correction 
ijRichards et al.ll2006 l for unobscur ed quasars. The total un- 
certainty introduced is ~ 0.5 dex ()Liu et al.ll2009l l. This is 
smaller than our smallest bin size (one dex) when we com- 
pare the scatter in the Mbh — o" relation with respect to 
nuclei luminosity. 

In order t o ob tain the nuclear lumino sity for the 
IGiiltekin et~aLl l|2009l) and I Graham et al] (|201ll ) samples, we 
select galaxies with nuclear luminosity measured in the soft 
X-Ray band by the Chandra X-ray observatorv llPellegrinil 
l2010l . l2005l : IZhang et al]|2009l : iGonzalez-Martm et al.ll2009^ 
At the lowest luminosity in our sample (~ 2 x 10^* erg s~^), 
the central X-ray source has luminosity compar able to the 

138-40 gj.g g-i) iKing et al.1 



in the NASA/IPAC Extraglactic Database (NEOfl Because 
these galaxies have nuclear luminosity(~ 2 x 10*^ erg s~^) 
much higher than the X-ray binaries, the X-ray luminos- 
ity of the galaxy is dominated by the AGN. Thus, we 
associate the X-ray luminosity of the galaxy with that 
of the AGN. Most of the X-ray data are obtained from 
XMM-Newton (iBianchi et al.1 l2009l : iMarkowitz et al.1 l2009l : 
iNandra et allbOOTl ), exce pt for Mrk202 an d IC 120, which 
are obtained frorti ASC A (|Ueda et al.l l2005') and BeppoSAX 
jVerrecchia et al1l2007l ) separately (detailed properties and 
references see Table 2). In order to obtain the X-ray lu- 
minosity from X-ray fluxes, we estimate distances assum- 
ing Ho = 73 km s-^Mpc-\ n„ = 0.27, = 0.73 for 
z > 0.01 galaxies, with redshift measurements obtained from 
NED, which compiled multiple consistent redshift measure- 
ments for each galaxy. For 2 < 0.01 galaxies, we obtain 
distances from the Extragalactic Distance Database, which 
gives updated bes t distances for galaxies within 3000 km s~^ 
(|Tu11v et al1l2009l ). 

To convert the X-ray luminosity to bolometric lumi- 
nosity, we first convert the X-ray luminosity of difi'erent 
bands in the literature to luminosity in the band 2 — 10 keV 
assuming an energy index of —1 {i^fi, = constant) with 
an uncertainty factor of 2 (Martin Elvis, private com- 
munication). Then, we convert the X-ray luminosity to 
bolometric luminosity by the bolometric correction for 
AGN s: Lb oi/Lx = 15.8, with an uncertainty of ~ 0.3 dex 
(Ho ,20091 ). Therefore, the nuclear bolometric luminosity is 
calculated as: 



Lboi ~ 15.8I/J5 



ln(10/2) 



(1) 



X-ray binary population(~ 10'^ _ 

l200lh . The use of Chandra data is therefore essential be- 
cause it has sufficient angular resolution to isolate galactic 
nuclear from bright X-ray binaries. For the reverberation- 
based sample, we select those with known X-ray luminosity 



ln(£;2/-Ei) ' 

where E2 and Ei represent the upper and lower bound of 
the observed X-ray band. 

We include the properties of the BH and host spheroids 
for our dynamically-based and reverberation-based sample 
in Table 1 and Table 2 separately. 

In reality the correction factor, Lboi/Lx, depends on 
the nuclea r luminos i ty: low luminosity AGNs tend to be "X- 
ray-loud" (|Hdl 19991 . [20091 ). In other words, the lower X-ray 
nuclear luminosity corresponds to an even lower bolometric 
nuclear luminosity and vice versa. Note that this additional 
complexity does not mix the order of galaxies with respect 
to their nuclei luminosity. As we compare galactic properties 
of lower nuclear luminosity galaxy to those of higher nuclear 
luminosity galaxies, without computing the exact bolomet- 
ric luminosity, our conclusion is not affected by assuming 
a constant bolometric correction factor. Because the bolo- 
metric cor rection fa ctor only introduces an uncertainty of 
- 0.3 dex (|Holl2009l ), the total uncertainty of the bolomet- 
ric luminosity is less than the smallest bin width, one dex. 
Therefore, we expect our results to be largely unaffected by 
the uncertainties in nuclear bolometric luminosity. 

In total, we have 38 galaxies with dynamically mea- 
sured BH mass and 17 galaxies with reverberation-mapping 
BH mass measurements in our sample. For the dynamically- 
based sample, the range of the nuclear luminosity is limited 
because the nuclear luminosity is difficult to measure for low 
nuclear luminosity galaxies and the BH mass is difficult to 

^ http://ned.ipac.caltech.edu/ 
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Table 1. Properties of Dynamically- based sample 



Name 




fj 


Co- 


log Mbh 


^log Mbh 


log[-^] 


ref 






km s~^ 


km s^^ 








(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


IC1459 (a) 


340 


17 


9.45 


0.34 


-5.52 


I'm \ 

(1) 


M31 (a 




160 


8 


8.18 


0.09 


-8.67 


(5) 


M32 (a 




75 


3 


6.49 


0.08 


-7.48 


(1) 


M60 (a 




385 


19 


9.32 


0.37 


-8.17 


(1) 


M81 (a 




143 


7 


7.90 


0.06 


-4.15 


(4) 


M84 (a 




296 


14 


9.18 


0.35 


-6.62 


(1) 


M87 (a 




375 


18 


9.56 


0.37 


-5.70 


(1) 


N524 (b) 


253 


25 


8.92 


0.10 


-7.27 


I'm \ 

(1) 


N1194 


c) 


148 


24 


7.82 


0.05 


-2.17 


(6) 


N2787 


a) 


189 


9 


7.63 


0.04 


-6.35 


(3) 


N2974 


b) 


227 


23 


8.23 


0.08 


-4.85 


(1) 


N3115 


a) 


230 


11 


8.98 


0.13 


-7.19 


(1) 


N3227 


a) 


133 


12 


7.18 


0.32 


-5.12 


(1) 


N3245 


a) 


205 


10 


8.34 


0.10 


-4.52 


(2) 


N3377 


a) 


145 


7 


8.04 


0.04 


-6.74 


(1) 


N3379 


a) 


206 


10 


8.08 


0.33 


-6.90 


(1) 


N3384 


a) 


143 


7 


7.26 


0.02 


-6.09 


(1) 


N3414 


b) 


237 


24 


8.40 


0.07 


-7.19 


(1) 


N3585 


a) 


213 


10 


8.53 


0.08 


-6.54 


(1) 


N3607 


a) 


229 


11 


8.08 


0.14 


-6.22 


(1) 


N3608 


a) 


182 


9 


8.32 


0.14 


-7.05 


I'm \ 

(1) 


N4026 


a) 


180 


9 


8.32 


0.08 


-6.85 


(1) 


N4261 


a) 


315 


15 


8.74 


0.09 


-4.58 


(1) 


N4459 


a) 


167 


8 


7.87 


0.39 


-6.41 


(1) 


N4552 


b) 


252 


25 


8.68 


0.07 


-6.37 


(1) 


N4596 


a) 


136 


6 


7.92 


0.13 


-6.40 


(2) 


N4621 


b) 


225 


23 


8.60 


0.07 


-6.62 


(1) 


N5128 


a) 


150 


7 


7.65 


0.13 


-3.54 


(1) 


N5813 


b) 


239 


24 


8.85 


0.07 


-5.24 


(1) 


N5846 


b) 


237 


9 


9.04 


0.04 


-6.70 


(1) 


N6251 


a) 


290 


14 


8.78 


0.18 


-2.36 


(2) 


N7457 


a) 


67 


3 


6.61 


0.25 


-5.64 


(1) 


UGC3789 (c) 


107 


12 


7.05 


0.05 


-0.82 


(6) 



Notes. Column 1: galaxy name and references on Mbhi c Column 2: stellar velocity dispersion (cr). Column 3: measurement error on 
cr. Column 4: black hole mass (Mbh)- Column 5: measurement error on logMBH(A^0). Column 6: X-ray luminosity references. 
References. 

A/ bh and a measuremen ts ( Column 1.) 

(a) iGultekin et~all |20o3) (b) ICraham et"all HoTJ) (c) ICreene et all l|2010l ') 

Nu clei X-r ay luminosity measurements (Column 7.) 

(DiPelleg rini (2010); (2) IConzalez-Martm et"all l|2009l ): (3') IPellegrinil {200^ : (4') IZhang et al.l l|2009f ): (51 1^ et al.l l|201lh : (6) 
ICreene et al.i yOlOj) 



estimate dynamically for high nuclear luminosity galaxies. 
The reverberation mapping measurements are normalized 
by setting a constant virial coefflcient so that the rever- 
beration mapping-based and dynamically-based Mbh — cr 
relations agree. The assumption of a constant virial coef- 
ficient could potentially introduce a larger scatter for the 
reverberation-based sample. In addition, the virial coeffi- 
cient may depend on the nuclear luminosity. If so, choosing 
a constant virial coefficient may introduce a BH mass un- 
certainty that depends on the nuclear luminosity. This could 
affect our result. In order to mitigate this potential bias, we 
do not rely heavily on comparing galaxies across our two 
samples. 

With Mbh, o" and nuclear luminosity in hand, we are 
now in a position to consider the scatter in the Mbh — cr re- 



lation with respect to nuclear luminosity. In order to obtain 
the scatter, we perform a linear fit to the Mbh — cr rela- 
tion by mi nimizing ■y'^ for the dynamically-based sample, 
as done bv lTremaine et al] ([2002). Our best-fit parameters 
are consistent with those in the literature: log (Mbh /Mq) = 
8.27 -I- 4.05 log((T/200 km s^^), with an intrinsic scatter of 
0.27. The scatter between the measured Mbh for an indi- 
vidual galaxy and the Mbh — cr relation, in units of the 
uncertainty, is calculated as: 



scatter = 



AMe 



(2) 



where AMbh is the difference between the measured BH 
mass and the BH mass corresponds to the fitted Mbh — cr 
relation, eM is the measurement uncertainty of BH mass. 
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Table 2. Properties of Reverberation-based sample 



Name a logMen eiosA/sH '"Strf^l '''^^ 

km s~i km s^^ Mq Mq 









yv 


V>) 


(P,\ 
V°) 




3C 120 


162 


20 


7.72 


0.23 


-0.18 


(5) 


Ark 120 


221 


17 


8.15 


0.11 


-0.60 


(4) 


Mrk 79 


130 


12 


7.70 


0.16 


-0.99 


(6) 


Mrk 110 


91 


7 


7.38 


0.14 


0.14 


(3) 


Mrk 202 


78 


3 


6.13 


0.22 


-0.58 


(2) 


Mrk 279 


197 


12 


7.52 


0.15 


-0.52 


(3) 


Mrk 590 


189 


6 


7.66 


0.12 


-0.46 


(7) 


Mrk 1310 


84 


5 


6.33 


0.17 


-0.30 


(1) 


Mrk 1383 


217 


15 


9.09 


0.16 


-1.21 


(3) 


N3227 


136 


4 


7.60 


0.24 


-2.06 


(5) 


N3516 


181 


5 


7.61 


0.18 


-1.13 


(5) 


N3783 


95 


10 


7.45 


0.13 


-0.10 


(3) 


N4051 


89 


3 


6.18 


0.19 


-1.05 


(8) 


N4151 


97 


3 


7.64 


0.11 


-1.46 


(4) 


N4253 


93 


32 


6.23 


0.30 


0.11 


(3) 


N4593 


135 


6 


6.97 


0.14 


-1.07 


(4) 


N5548 


195 


13 


7.80 


0.10 


-0.99 


(4) 


N6814 


95 


3 


7.25 


0.12 


-1.52 


(4) 


N7469 


131 


5 


7.06 


0.11 


-0.32 


(9) 



Notes. Column 1: galaxy name. Column 2: stellar velocity dispersion (cr). Column 3: measurement error on a. Column 4: black hole 

mass (J\/bh). Column 5: measurement error on logMBuC^©)- Column 6: X-ray luminosity references. 

References. 

Nu clei X-ray lumino sity measurements and instruments (Column 7.) 

(DlVerrecchia et al. | ( 2007). BeppoSAX; {2a) Bianchi et al. (2009), XMM; ('2b') iMarkowitz et al.l (|2009|) . XMM; (2c') lNandra et al.l 
l|2007l '). XMM: f3') IUeda et al.l (l2005h . ASCA 

Local galaxy (z<0.01) distance measurements: N3227, N3516, N3783, N4051, N4151, N4593: lTuilv et al.l ll2009l) : 

ecr is the measurement uncertainty of the velocity dispersion 
and b is the slope of the A/bh — o" relation. For a fixed slope 
of the Mbh — o" relation, the scatter in BH mass determines 
the scatter in a [scatter bh — b x scatter^). Measuring the 
scatter in BH mass makes it more convenient to compare 
with the liter ature, where intrinsic scatter in BH mass is 
discussed (eg. iTremaine et al]|2002l ). 



3 RESULTS 

In Figure 1 we plot the Mbh — cr relation for our sample. In 
order to compare the scatter in the Mbh — o" relation with 
the nuclear luminosity, we color code the symbols by the 
nuclear luminosity level, log Ltoi/LEdd- Because the rever- 
beration mapping method has a higher uncertainty in BH 
mass (because a constant virial coefficient is adopted), the 
scatter in the reverberation-based sample is slightly larger 
than that for the dynamically-based sample. It is apparent 
already from this figure that there is no strong correlation 
between Mbh — cr scatter and nuclear luminosity within each 
sample. 

In order to investigate the dependence of the scatter on 
nuclear luminosity quantitatively, we show in Figure 2 the 
scatter versus the nuclear luminosity (top panel) and the 
average scatter in bins of the nuclear luminosity level (bot- 
tom panel). To keep track of the standard deviation of the 
scatter at each nuclei luminosity bin, we plot the standard 
deviation in each bin as vertical dash error bars. We choose 



log[Lboi/LEdd] 




-0.4 -0.2 0.2 0.4 

log[cr/(200 km s-i)] 

Figure 1. A/bh — o" relation for galaxies with classical bulges and 
dynamically- based BH masses (stars), and reverberation-based 
masses (circles). The color bar indicates the nuclear luminosity 
levels, log Li,oi/LEdd- This figure demonstrates visually that there 
is no strong correlation between the scatter and the nuclear lu- 
minosity. 
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Figure 2. Scatter (top panel) and binned scatter (bottom panel) 
in the Mbh — o" relation as a function of galaxy nuclear luminosity. 
Measurement error is taken to be the sum of the measurement un- 
certainty in Mbh ^^nd a, where the measurement uncertainty in a 
is scaled to Mbh (see scatter definition in equation (2)). Galaxies 
with dynamical BIf measurements are represented by stars, and 
those with reverberation measurements are represented by cir- 
cles. The solid vertical error bars indicate the error on the mean 
scatter in each bin; the dashed vertical error bars indicate the 
standard deviation of the scatters in each bin; the horizontal er- 
ror bars indicate the standard deviation of the nuclear luminosity 
levels in each bin. It is clear from this figure that the scatter is 
independent of nuclear luminosity. 



This result can be seen another way by considering the 
timescales that govern the growth of the BH mass and the 
dynamical time. We define the BH growth timescale as: 

Mbh 



Mb 



= 4 X 10^ 



e 



yr, 



(3) 



where e is the radiative efficiency, which we take to be 0.1, 



= 3.5 X 10*( 



i4/-A£bh ■ 







M, 



Lq is the Eddington luminosity. 



The dynamical time is defined via: 



a 



(4) 



where Re is the effective radius of the galaxy. Because we 
are interested in the two timescales when the BH is accreting 
rapidly, we set an arbitrary threshold (4 x lO^^Lfidd) and 
select galaxies with nuclear bolometric luminosity higher 
than this value. Our result is independent of this thresh- 
old. This leaves us with 3 dynamically-based measurements, 
and all the reverberation-based measurements. Then, to cal- 
culate tdyn, we obtain the effective ra.dii of these galax- 
ies from the literature (ISani et al.ll201ll: iBentz et aLl l2009l : 



iLauer et al]|2007l : iMarconi fc Hunt! 120031 ). For the galaxies 
with no measured effective radii in the literature (N1194, 
UGC3789, Mrk 202 and N4253), we use the galaxy radius 
calculated as the following. We estimate the radii of the 
galaxy by multiplying the angular radii by the angular-size 
distances. We obtain the angular radii from the 2MASS 
isophotal measurements with reference level of the radii 
set at 20 K — band magnitude arcsec"^. The angular-size 
distance is calculated assuming Ho = 73 km s"^Mpc"^ 
f^m ~ 0.27, Oa = 0.73 as the four galaxies are all z > 0.01 
galaxies. The redshift of the galaxies are obtained from 
NED, which compiled multiple consistent redshift measure- 
ments for each galaxy. 

In Figure 3 we plot the ratio of tacc to tdyn as a function 
of the nuclear luminosity in units of the Eddington luminos- 
ity. We fit a straight line and find the slope of the line is 
—0.94 and the ratio approaches unity when the nuclear lu- 
minosity approaches the Eddington limit. As discussed in 
section 2, the bolometric correction factor depends on the 
nuclei luminosity. By fixing a constant correction factor, the 
nuclear bolometric luminosity is probably underestimated 
at high luminosity by roughly factor of 2. This causes tacc 
to be overestimated at high luminosity by the same factor. 
Thus, plotting the "true" tacc versus the "true" luminosity, 
which are overestimated and underestimated separately by 
the same factor, the points in Figure 3 should shift to the 
left along our observed trend. Therefore, the trend in Figure 
3 is robust against the uncertainties in nuclear luminosity. 



bins so that the number of galaxies per bin is comparable, 
in order to minimize Poisson noise. Specifically, the number 
of galaxies in each bin are 13, 13, 12, 9, 8. Notice that we 
do not mix the dynamically-based and reverberation-based 
samples. The scatter and the standard deviation remain ap- 
proximately constant as the nuclear luminosity increases. 
This is the basic result of this article. 



4 DISCUSSION 

The relationship between scatter in the A4"bh — cr relation and 
the accretion rate of the BH puts interesting constraints on 
how galaxies evolve in the Mbh — cr plane. For example, if 
BH growth and host stellar spheroid growth are uncorrelated 
then we would expect that galaxies with rapidly accreting 
BHs to lie systematically off of the A^bh — cr relation. In con- 
trast, if BHs and spheroid coevolve, then we would expect 
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-2 -1.5 -1 -0.5 

log[Lboi/LEdd] 



Figure 3. Ratio of the instantaneous accretion timescale to dy- 
namical timescale as a function of bolometric nuclear luminosity 
in units of the Eddington luminosity. The black line indicates 
the linear fit of the relation, which has a slope of -0.94. The fact 
that the ratio of timescales approaches unity as the luminosity 
approaches Eddington implies that BH mass and a coevolve. 



no correlation between scatter in Mbh — o" and BH activity, 
which is precisely what we observe. Our results therefore fa- 
vor the scenario wherein spheroid and BHs coevolve along 
the Mbh — o" relation. 

In support of our conclusion, we also find that the dy- 
namical time of the galaxy becomes comparable to the BH 
accretion timescale when the nuclear luminosity approaches 
the Eddington luminosity (Figure 3). We can understand 
how galaxies populate the space in Figure 3 in light of our 
results. Galaxies cannot be in the upper right region of Fig- 
ure 3 because otherwise the BH would grow much slower 
than a, resulting in a correlation between scatter and nuclear 
luminosity, which is not observed. We can also understand 
why there are no galaxies in the lower left region of Figure 
3: there is no limit to how long tacc can be (because the 
accretion rate can be arbitrarily close to zero), while there 
is a limit to how small tdyn can be. So tacc/tdyn can increase 
at lower L/LEdd- For systems experiencing rapid growth in 
the BHs, the growth timescale of the BH is comparable to 
the dynamical time, supporting the idea that galaxies evolve 
along the Mbh — o" relation. 

As a further extention, we also consider the intrinsic 
scatter of galaxies with Mbh measu red by the virial method 
IXiao et al.ll201ll : IShen et al.ll201ll ) for high redshift AGN. 
We find that the intrinsic scatter of those is similar to the 
intrinsic scatter of the galaxies with BH masses measured by 
the reverberation-mapping method. The virial method uses 
the radius-luminosity relation to estimate the radius of the 
broad-line region, and then uses the reverberation mapping 
method formalism to measure BH masses. Thus, the mass 
of the BH has an even bigger uncertainty. Barring these 
additional uncertainties, this suggests that high-z galaxies 
may also evolve along M — a relation. 

Our result is inconsistent with models that predict a 



non-causal origin of BH mass-galaxy property relations. As 
discussed in ljahnke fc Macciol ( 201 ih . in a non-causal origin 
model, the BH mass growth rate is uncorrelated with the 
growth of a, and the BH mass-galaxy property relation con- 
verge only through the central limit theorem. Thus, when a 
given BH is growing its mass at high nuclear luminosity, a 
does not "catch up" until after several merger events. Such 
a model would therefore predict the scatter in Mbh — cr to 
be larger for higher nuclear luminosity, in contradiction to 
our results. 

On the other hand, our conclusion is consistent with 
the scenario emerging from simulations that episodes of 
major spheroid growth and BH growth occur on sim- 
ilar timescales via mergers. Simulations find that the 
epoch of rapid BH accretion is limited by AGN feed- 
back to be ~ 100 Myr, which is similar to the dynami- 

) — II 1 I 1 

cal time (eg. KaufFmann & Haehnelt 2000; Di Matteo et alj 

2005; Sijacki ct al. 2007; Hopkins et al. 2009; Hopkins 201il; 
iBlecha et al.ll201ll ). This suggests that the BH mass and a 
change concurrently. Mo dels that appeal to fue ling of the 
BH by recycled gas (e.g.. ICiotti fc Ostrikeij[2007[ ) must sat- 
isfy our constraint that galaxies evolve along Mbh — o" even 
during high accretion rate phases. 

In addition, our result reinforces the assumptions in var- 
ious studies. For instance, it shows that BH grows simultane- 
ously with t he potential well form ation in mergers, which is 
assumed bv lShankar et all (|2009l ) to investigate the cosmo- 
logical evolution of the Mbh — o" relation. The self-regulated 
growth of supermassive BH is assumed when investigating 
the link betwe en quasars and the red galaxy population by 
iHopkins et al.l (2006), an d when constraining th e accretion 
history of massive BHs bv lVolonteri et all (|2006l ). Moreover, 
it constrains the total tim escale of the episod ic random ac- 
cretion model proposed by I Wang et al. I (|2009t ) to be similar 

to tdyn- 

In our dynamically-based sample we only have one 
galaxy with nuclear bolometric luminosity higher than 
O.lLEdd, and three galaxies with nuclear bolometric lumi- 
nosity higher than O.OOlI/sdd. We can expect more strin- 
gent constraints on the coevolution of galaxies and BHs as 
additional maser-based BH masses are obtained for higher 
nuclear luminosity galaxies. 



ACKNOWLEDGMENTS 

We thank Martin Elvis, Pepi Fabbiano, Silvia Pellegrini, 
Junfeng Wang, Yue Shen, Laura Blecha and Phil Hopkins for 
helpful conversations. This work was supported in part by 
NSF grant AST-0907890 and NASA grants NNX08AL43G 
and NNA09DB30A. 



REFERENCES 

Beifiori A., Courteau S., Corsini E. M., Zhu Y., 2011, ArXiv 
e-prints 

Bentz M. C., Peterson B. M., Netzer H., Pogge R. W., 

Vestergaard M., 2009, ApJ, 697, 160 
Bianchi S., Guainazzi M., Matt G., Fonseca Bonilla N., 

Ponti G., 2009, AAP, 495, 421 



© 0000 RAS, MNRAS 000, 000-000 



Coevolution of Supermassive Black Holes and their Host Spheroids 7 



Blecha L., Cox T. J., Loeb A., Hernquist L., 2011, MNRAS, 
412, 2154 

Burkert A., Silk J., 2001, ApJL, 554, L151 
Ciotti L., Ostriker J. P., 2007, ApJ, 665, 1038 
Di Matteo T., Springcl V., Hernquist L., 2005, Nature, 433, 
604 

Fabian A. C, 1999, MNRAS, 308, L39 

Ferrarese L., Merritt D., 2000, ApJL, 539, L9 

Gebhardt K., Bender R., Bower G., Dressier A., Faber 
S. M., Filippenko A. V., Green R., Grillmair C., Ho L. C., 
Kormendy J., Lauer T. R., Magorrian J., Pinkney J., 
Richstone D., Tremaine S., 2000, ApJL, 539, L13 

Gonzalez-Martin O., Masegosa J., Marquez L, Guainazzi 
M., 2009, ApJ, 704, 1570 

Graham A. W., Onkcn C. A., Athanassoula E., Combes F., 
2011, MNRAS, 412, 2211 

Greene J. E., Peng C. Y., Kim M., Kuo C.-Y., Braatz J. A., 
Violette Impellizzeri C. M., Condon J. J., Lo K. Y., Henkel 
C, Reid M. J., 2010, ApJ, 721, 26 

Giiltekin K., Richstone D. O., Gebhardt K., Lauer T. R., 
Tremaine S., Aller M. C, Bender R., Dressier A., Faber 
S. M., Filippenko A. V., Green R., Ho L. C, Kormendy 
J., Magorrian J., Pinkney J., Siopis C, 2009, ApJ, 698, 
198 

Giiltekin K., Tremaine S., Loeb A., Richstone D. O., 2011, 
ApJ, 738, 17 

Hirschmann M., Khochfar S., Burkert A., Naab T., Genel 

S., Somerville R. S., 2010, MNRAS, 407, 1016 
Ho L. C, 1999, ApJ, 516, 672 
— , 2009, ApJ, 699, 626 
Hopkins P. F., 2011, ArXiv e-prints 

Hopkins P. F., Hernquist L., Cox T. J., Keres D., Wuyts 

S., 2009, ApJ, 691, 1424 
Hopkins P. F., Hernquist L., Cox T. J., Robertson B., 

Springel V., 2006, ApJS, 163, 50 
Jahnke K., Maccio A. V., 2011, ApJ, 734, 92 
Johansson P. H., Naab T., Burkert A., 2009, ApJ, 690, 802 
KaufTmann G., Haehnelt M., 2000, MNRAS, 311, 576 
Kauffmann G., Heckman T. M., Tremonti C, Brinchmann 

J., Chariot S., White S. D. M., Ridgway S. E., Brinkmann 

J., Fukugita M., Hall P. B., Ivezic Z., Richards G. T., 

Schneider D. P., 2003, MNRAS, 346, 1055 
King A. R., Davies M. B., Ward M. J., Fabbiano G., Elvis 

M., 2001, ApJL, 552, L109 
Kormendy J., Bender R., Cornell M. E., 2011, Nature, 469, 

374 

Kormendy J., Gebhardt K., 2001, in American Institute of 
Physics Conference Series, Vol. 586, 20th Texas Sympo- 
sium on relativistic astrophysics, J. C. Wheeler & H. Mar- 
tel, ed., pp. 363-381 
Kormendy J., Richstone D., 1995, ARA&A, 33, 581 
Lauer T. R., Gebhardt K., Faber S. M., Richstone D., 
Tremaine S., Kormendy J., Aller M. C, Bender R., 
Dressier A., Filippenko A. V., Green R., Ho L. C, 2007, 
ApJ, 664, 226 

Li Z., Garcia M. R., Forman W. R., Jones C, Kraft R. P., 
Lai D. v., Murray S. S., Wang Q. D., 2011, ApJL, 728, 
LIO-I- 

Liu X., Zakamska N. L., Greene J. E., Strauss M. A., Krolik 

J. H., Heckman T. M., 2009, ApJ, 702, 1098 
Marconi A., Hunt L. K., 2003, ApJL, 589, L21 
Markowitz A., Reeves J. N., George I. M., Braito V., Smith 



R., Vaughan S., Arevalo P., Tombesi F., 2009, ApJ, 691, 

922 

Murray N., Quataert E., Thompson T. A., 2005, ApJ, 618, 
569 

Nandra K., O'Neill P. M., George I. M., Reeves J. N., 2007, 

MNRAS, 382, 194 
Pellegrini S., 2005, ApJ, 624, 155 
— , 2010, ApJ, 717, 640 

Peng C. Y., 2007, ApJ, 671, 1098 

Reyes R., Zakamska N. L., Strauss M. A., Green J., Krolik 
J. H., Shcn Y., Richards G. T., Anderson S. F., Schneider 
D. P., 2008, AJ, 136, 2373 

Richards G. T., Strauss M. A., Fan X., Hall P. B., Jester 
S., Schneider D. P., Vanden Berk D. E., Stoughton C, 
Anderson S. F., Brunner R. J., Gray J., Gunn J. E., Ivezic 
Z., Kirkland M. K., Knapp G. R., Loveday J., Meiksin 
A., Pope A., Szalay A. S., Thakar A. R., Yanny B., York 
D. G., Barentine J. C, Brcwiugton H. J., Brinkmann J., 
Fukugita M., Harvanek M., Kent S. M., Klcinman S. J., 
Krzesihski J., Long D. C, Lupton R. H., Nash T., Neilsen 
Jr. E. H., Nitta A., Schlegel D. J., Snedden S. A., 2006, 
AJ, 131, 2766 

Richstone D., Ajhar E. A., Bender R., Bower G., Dressier 
A., Faber S. M., Filippenko A. V., Gebhardt K., Green 
R., Ho L. C, Kormendy J., Lauer T. R., Magorrian J., 
Tremaine S., 1998, Nature, 395, A14-I- 

Robertson B., Cox T. J., Hernquist L., Franx M., Hopkins 
P. F., Martini P., Springel V., 2006, ApJ, 641, 21 

Sani E., Marconi A., Hunt L. K., Risaliti G., 2011, MNRAS, 
413, 1479 

Shankar F., Bernardi M., Haiman Z., 2009, ApJ, 694, 867 
Shen Y., Richards G. T., Strauss M. A., Hall P. B., 

Schneider D. P., Snedden S., Bizyaev D., Brewington H., 

Malanushenko V., Malanushenko E., Oravetz D., Pan K., 

Simmons A., 2011, ApJS, 194, 45 
Sijacki D., Springel V., Di Matteo T., Hernquist L., 2007, 

MNRAS, 380, 877 
Silk J., Nusser A., 2010, ApJ, 725, 556 
Silk J., Rees M. J., 1998, AAP, 331, LI 
Tremaine S., Gebhardt K., Bender R., Bower G., Dressier 

A., Faber S. M., Filippenko A. V., Green R., Grillmair 

C, Ho L. C, Kormendy J., Lauer T. R., Magorrian J., 
Pinkney J., Richstone D., 2002, ApJ, 574, 740 

Tully R. B., Rizzi L., Shaya E. J., Courtois H. M., Makarov 

D. I., Jacobs B. A., 2009, AJ, 138, 323 

Ueda Y., Ishisaki Y., Takahashi T., Makishima K., Ohashi 

T., 2005, ApJS, 161, 185 
Verrecchia F., in't Zand J. J. M., Giommi P., Santolamazza 

P., Granata S., Schuurmans J. J., Antonelli L. A., 2007, 

AAP, 472, 705 

Volonteri M., Salvaterra R., Haardt F., 2006, MNRAS, 373, 
121 

Wang J.-M., Hu C, Li Y.-R., Chen Y.-M., King A. R., 
Marconi A., Ho L. C, Yan C.-S., Staubert R., Zhang S., 
2009, ApJL, 697, L141 

Woo J.-H., Treu T., Barth A. J., Wright S. A., Walsh J. L., 
Bentz M. C, Martini P., Bennert V. N., Canalizo G., Fil- 
ippenko A. v., Gates E., Greene J., Li W., Malkan M. A., 
Stern D., Minezaki T., 2010, ApJ, 716, 269 

Wyithe J. S. B., Loeb A., 2003, ApJ, 595, 614 

Xiao T., Barth A. J., Greene J. E., Ho L. C, Bentz M. C, 
Ludwig R. R., Jiang Y., 2011, ArXiv e-prints 



© 0000 RAS, MNRAS 000, 000-000 



8 Li, Conroy and Loeb 



Zakamska N. L., Strauss M. A., Krolik J. H., Collinge M. J., 
Hall P. B., Hao L., Heckman T. M., Ivezic Z., Richards 
G. T., Schlegel D. J., Schneider D. P., Strateva I., Vanden 
Berk D. E., Anderson S. F., Brinkmann J., 2003, AJ, 126, 
2125 

Zhang W. M., Soria R., Zhang S. N., Swartz D. A., Liu 
J. F., 2009, ApJ, 699, 281 



© 0000 RAS, MNRAS 000, 000-000 



